Collisional energy losses of heavy quarks in relativistic nuclear collisions at energies 
available at the BNL Relativistic Heavy Ion Collider 
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We investigate the energy losses of b and c quarks at hard coUisions with gluons and light quarks 
in quark-gluon plasma produced in central Au + Au collisions at RHIC energy (-^sjvjv = 200 GeV). 
As the couphng constant at hard collision in the perturbative QCD is limited by Landau pole at 
momenta = A^, we use the analytic models of QCD, where the Landau pole is absent. We 
calculate the energy losses of 6 and c quarks in analytic model of QCD at > A^. At calculations 
we use the effective quasiparticle model. We use at calculations of total nuclear modification factor 
R''j^^ the fragmentation of b and c quarks into heavy hadrons with accounting theirs energy losses. 
We show that this factor gives the considerably suppression of heavy B and D mesons in the region 
of middle meanings of p±. The experimental data 01 show, that nonphotonic electrons have the 
analogous suppression in the same region of p±. 

PACS numbers; 12.38.Mh, 24.85.+p, 25.75.-q 



I. INTRODUCTION 

Heavy quarks, unlikely light quarks do not come to an 
equilibrium with the surrounding matter and may there- 
fore play an important role in the search for properties 
of this matter. Produced in hard collisions, their initial 
momentum can be inferred from pp collisions. The de- 
viation of the measured heavy mesons p± distribution 
in AA collisions from that measured in pp collisions is 
usually quantified by the nuclear modification factor 



Raa = 



dNAu 



+Au 



Taa dop+p 



(1) 



where dNAu+Au is the differential electrons yield from 
heavy meson decays in Au + Au collisions and dap+p is 
corresponding cross section in p+p collisions in any given 
p± bin. 

In the RHIC experiments the heavy mesons have not 
yet been directly measured. Both the PHENIX Q 
and STAR 01 collaboration have observed only single 
nonphotonic electrons, which have been created in the 
semileptonic decay of heavy mesons. Thus experimen- 
tally it is not easy separate between charm and bottom 
hadrons. The QCD calculations in the fixed order -|- 
next to leading logarithm (FONLL) predict ^ the theo- 
retical uncertainty for the charm and bottom quarks p± 
distribution in pp-collisions at ^/snn = 200 GeV. Ap- 
parently, that above p± = 4 GeV/c the electrons from 
bottom mesons dominate the spectrum. However, un- 
certainty here is considerable. The nonphotonic electron 
yield at RHIC measurement exhibits unexpectedly large 
suppression in central Au + Au collisions at high p± . The 
observed values Raa ~ 0.2 are smaller than originally 
expected. The theoretical approaches based on pertur- 
bative QCD (pQCD) calculation give much larger values 
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and it has been doubt, whether pQCD is the right de- 
scription of this interaction. One should note that for 
light quarks the radiative energy loss is the most im- 
portant. The heavy quarks have a smaller amount of 
radiative energy loss due to their large mass [1] , . 

The asymptotic freedom in QCD enable one to ap- 
ply the standard perturbative formulas at large momenta 
transferred. However a number of phenomena are be- 
yond such calculation. Moreover, the perturbative the- 
ory contain unphysical singularities in the expression for 
the running coupling, for example the ghost poles and 
also the unphysical cutes. Therefore by number of au- 
thors were considered the analytic models of invariant 
charge in QCD. The condition of analyticity follows from 
the general principles of local Quantum Field Theory. 
The analytic models contain no unphysical singularities 
and running coupling may be determined in the spacelike 
(Euclidean) region and in the timelike region by analytic 
continuation from spacelike region. In this work we in- 
vestigate the collisional energy losses of b and c quark in 
quark-gluon plasma (QGP) produced in central Au + Au 
collisions at RHIC energy. We use here the formulas for 
analytic invariant charge in QCD in the spacelike region. 
We show that the nuclear modification factor R'a'a ob- 
served in the experiments is determined by hard collisions 
with QCD analytic charge aan{Q^) which have the Lan- 
dau pole at = A^ in the perturbative QCD, but have 
no of this pole in analytic QCD. The vicinity = A^ one 
should consider as boundary of hard collisions in QCD. 

In Sec.|TT]we give the brief description of analytic mod- 
els results, which we use in this paper, and corresponding 
reference. 

In Sec. mil we give the brief description of effective 
quasiparticle model and physical characteristic of ex- 
panding plasma. We calculate here energy losses AE/E 
of b and c quarks at collisions with gluons and light 
quarks in plasma, using the analytic coupling in QCD. 
We calculate also the summary losses of b and c quarks 
for various energies while they transverse the expanding 
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plasma. 

In Sec. IIVI we calculate the common nuclear modifi- 
cation factor R''j^^ for b and c quarks, which fragmen- 
tized into B and D mesons. We use the fragmentation 
function 6 and c quarks with accounting of theirs energy 
losses. We take into account the relative B and D con- 
tribution into electrons from theirs semileptonic decays. 
We compare this modification factor with experimental 
data. 

In Sec.|V]- Conclusion. 



(AIC) at the higher loop levels was also investigated in 
the works 3, Q on the basic of the Kallen - Lehmann 
type integral representation. It was shown that AIC pos- 
sesses the higher loop stability, i.e. the curves corre- 
sponding to the two and three-loop levels are almost in- 
distinguishable. In the present work we use analytic in- 
variant charge in the form Q for investigation of energy 
losses b and c quarks at collisions with gluons and light 
quarks in the quark-gluon plasma produced in central 
Au + Au collisions at RHIC energy. 



II. ANALYTIC MODELS FOR THE QCD 
RUNNING COUPLING 

The effective coupling constant a^iQ^) is given in the 
lowest approximation (one loop) by the famous asymp- 
totic formulas 



(11 



2jV/ 
3 



)ln. 



(2) 



where A is dimensional constant of QCD. This coupling 
constant have divergence at — A^. That is related 
to an infrared Landau pole in the running coupling, that 
is nonperturbative issue. In the two-loop case we have 
besides pole also the unphysical cut owing to the fac- 
tor of type In(ln^). In consequence of these difficul- 
ties the number of authors use the analytic models of 
invariant charge in QCD. For example, some of them; 
[e*!, [7], ¥], 9]. The Kallen - Lehman spectral represen- 
tation satisfied to condition of analyticity and it applied 
to the " analytization" of the renormalization group (RG) 
equations. We describe here briefly the analytic model 
proposed in the works , @ • Here the solution is derived 
as the solution of the RG equation with the (3 function 
"improved" by the analytization procedure. This corre- 
spond to demand of analyticity of right hand side of RG 
equations as a whole before their solution by use of the 
Kallen - Lehmann integral representation. At the Z-loop 
level the analytic running coupling ain(g^) is defined as 
the solution of equation: 



dlnai'2(g^) 
din 



The value a"an{q^) at the one- loop level has here the 
form: 



iYBMHew^'}. 



J=0 



(3) 



a 



(1) 



47r(z- 1) 



(4) 



where z — -jj, = 11 — This invariant charge 
have no of the Landau pole in the vicinity of the point 
z = 1, namely ain(z = 1 4- e) ~ f^(l ~ f) ^.nd it de- 
scribe also the asymptotic for large z. The infrared en- 
hancement at small z is in agreement with the Schwinger 
- Dyson equations fld\. The analytic invariant charge 



III. COLLISIONAL ENERGY LOSSES OF b AND 
c QUARKS IN QUARK-GLUON PLASMA 

By analogy with the works we use the effec- 

tive quasiparticle model for calculation of physical char- 
acteristics of expanding plasma. We assume the hot 
glue production at the early stage. This is caused by 
the relatively large gg cross section in comparison with 
qg and qq ones [l3| . In the work were calculated 
the initial temperature Tq and time tq at y^J/vw = 200 
GeV and also t(T) for evolution of the plasma phase 
for T — )■ Tc from above. For these estimations we used 
the formulas of the quasiparticle model for gluons and 
quarks densities Ug, n^, rig at temperature T at expand- 
ing of plasma [l3| ■ These densities depend from thermal 
masses mg{T) and mq[T) of gluons and quarks. These 
thermal masses are connected with the thermodynamic 
functions of medium and they are determined by effec- 
tive coupling (T) . They describe well the lattice data 
aX Tc < T < ATc. We use here the parametrization of 
effective two-loop coupling G^{T,Tc) from Ref. [l3]- We 
investigate in this work the central Au + Au collisions 
for nuclei modeled as a uniform cylinder with sharp edge 
and effective radius Reff — 6 fm, which correspond to 
T ~ 6 fm and T ~ 183 MeV This value of T is 

close to Tc = 177 MeV [11]. The calculation of the lep- 
ton spectrum includes initial b and c quarks distribution 
from perturbative QCD, heavy quarks Q energy losses, 
their fragmentation into heavy mesons Hq and electron 
decay spectrum of heavy mesons. This cross section is 
schematically written as [1] , [H] : 



E(fa{e 



E,(Pa{Q) 
dpi 



g)P(i?.^i?/)(g) 



(5) 



where the symbol denotes a generic convolution. The 
term /{Hq — >■ e) includes the branching ratio to elec- 
trons. The initial b and c p± distributions in pp collisions 
with account of theoretical uncertainty were computed at 
FONLL level in Ref. In this work we use these dif- 
ferential cross section ^y^p^ (look at Fig.[T]). 

The invariant cross section for collisions of heavy quark 
Q with momentum p and massless partons with momen- 
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turn k in plasma can be written in the form: 



da = dp'dk' 



4a2 \Ml 
4J 



AEp-k' 



(5(p + k - p' - k') X 

-Ifc'l) 



= dp'dk'F(p,k,p',k',mQ), 



(6) 



where p' and k' are the momenta of outgoing heavy 
quarks Q and hght partons in plasma, |Af p is the ma- 
trix element squared. The value J = y/ (PK)"^, where 
PK = Ep\k\ - pk = k{Ep - pcos^) and * the an- 

S—Al^ 

gle between vectors p and k. Thus J = — where 
S = + 2\k\{Ep - poos*). Using the formula dH) for 
cross section da, one should determine the relative energy 
loss per unit length (for Qg collisions) from integral: 



(ikn„(k) w 
— -— — ^ da, 

(27r)3 p 



(7) 



where uj = Ep — Ep/ (in t channel) and '^g(k) is the gluon 
density in the quasiparticle model. We calculate in these 
equations at first the fc'-integral, using three-dimensional 
6 function, and this gives k' = |k-|-p — p'|. For the second 
S function it is convenient to choose for massless partons 
the z axis along the direction p of the heavy quarks, and 
k contained in the yz plane [12] , [13 ■ We use the second 
5 function for calculation of the integral / over azimuthal 
angle (j) of the form: 

1 = J 5(^Ep~Ep, +k- y/XTWcospj . (8) 



The result of calculation is: 



- 



(9) 



In values A and B we change the variables from p' 
and cos 9 to variables t and w with corresponding Jaco- 
bian, and t = 2M^ - 2{EpEp' — pp' cos 6) . The formula 
of type © was used in Ref. fT?! for investigation of col- 
lisional energy loss of muons with mass M in the hot 
QED plasma. The value — has here the form: 
B^ — = —d}ijj^ -^hu) + c. We use in this paper the 
coefficients a, b, c from Ref. for the masses Mq of 
heavy quarks. The value B^ — is positive in interval 



< U! < CO- 



nant D 



max: 

6- 



where lo™^^ 



is 



0: 



b±VD 
2a^ 



The discrimi- 



D = -t[st + (s - MA)2] 



2 2i f 4fcsin\[' 



V P 



The condition D > gives < \t\ < 



{s-Mlf 



(10) 



How- 



ever we use here the analytic model of QCD with lower 
boundary of — t — in the vicinity of Landau pole t ~ 
where the QCD running coupling for hard collisions con- 
tain no unphysical singularity. We can evaluate now the 



relative energy loss for h and c quarks using the equa- 
tions ^ and ([7]) and the lj integral: 



jduj 



',\t\p[Ep{s-M^Q)-k{s + Ml)] 



- 



(S-M2)3 



We have thus the formula (for Qg collisions): 



(11) 



'Ed^ 



dk rtg (k) 
(27r)3 

(s -I- M2) 



d\t\al^{t)Ti\t\ X 



|t|n,.„ 

E(s-M'^) 



\M, 



Qa\ 



{s^APfE 



4J 



(12) 



where 4 J = 2(s - M2), and M is Mb or Mc - the mass of 
6 or c quark. The square of the matrix element in t chan- 



nel have the form llE 



-\t\ 



where s — s — AP 



In quasiparticle model the thermal distribution of glu- 
ons is ng{k) = , where k — yT, and 



cxp (V + 



^ (T) 
?2— )-l 



the thermal mass rUg (T) is expressed across the effective 
two - loop coupling G^{T, Tc) from Ref. [ij. From equa- 
tion (fT2| we have the condition (E — ^ ) > k. It 
is convenient to introduce the variable = E/p — cos^, 
i.e. the condition is E — ¥— > k. We have also 
s — 2yT^ip. From this condition and also from the con- 



dition ^ > \t[ 



we have k > 



1 



2M 



and 



6 > 



1 



Ep \ ' 2M 

in formula (1121) is 



The upper limit for variable \t\ 



1 



(13) 



We use for calculations the mass of heavy quarks Mf, = 
4.75 GeV, Mc = 1.3 GeV. For example, at p = 8 GeV we 
have Eb = 9.304 GeV and E^ ^ %A GeV. We have also 
Cimax = {E/p+l). We write down now the formula (fT^ 
for heavy quarks with the mass M with accounting of 
analytic invariant charge in the form Q: 



'Ed^ 



16TT2 /47r 
~A^pE \% 



E/p+l 



db_ 
26 



dyy- 



f (l 



(2!/5l)^ 



dx 

z 



z- 1 
z In z 



(14) 
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Table I: The relative energy losses for 6 and c quark at bg and 
eg collisions 



r(MeV) 


T(m-^) 






410.6 


0.405 


0.143 


0.49 


390 


0.474 


0.128 


0.433 


370 


0.543 


0.113 


0.379 


350 


0.613 


0.0993 


0.329 


330 


0.691 


0.086 


0.283 


310 


0.825 


0.0738 


0.239 


290 


0.98 


0.0618 


0.199 


270 


1.07 


0.0508 


0.162 


250 


1.19 


0.0404 


0.127 


240 


1.26 


0.0355 


0.115 


233 


1.31 


0.0324 


0.1 


220 


1.62 


0.026 


0.0806 


210 


1.92 


0.0211 


0.0656 


200 


2.36 


0.0163 


0.0505 


190 


3.08 


0.0112 


0.0347 


183 


4.13 


0.0073 


0.0226 



Table II: The relative energy losses for b and c quarks at bq 
and cq collisions 



r (MeV) r jm-') (|), 



410.6 


0.405 


0.0295 


0.1013 


390 


0.474 


0.036 


0.1222 


370 


0.543 


0.0402 


0.1351 


350 


0.613 


0.0454 


0.1511 


330 


0.691 


0.0502 


0.165 


310 


0.825 


0.055 


0.1803 


290 


0.98 


0.0612 


0.198 


270 


1.07 


0.051 


0.1627 


250 


1.19 


0.042 


0.133 


240 


1.26 


0.0381 


0.12 


233 


1.31 


0.0348 


0.1085 


220 


1.62 


0.029 


0.0903 


210 


1.92 


0.025 


0.0776 


200 


2.36 


0.0205 


0.0635 


190 


3.08 


0.0156 


0.0484 


183 


4.13 


0.0113 


0.035 



where ^imin 



= ^fl 



Bp V-^ ' 2M ^min — 25ipT 

here the analytic approach in QCD for investigate of en- 
ergy losses of heavy quarks at hard coUisions vifith light 
partons in plasma at temperature Tq > T > Tc- We in- 
vestigate the energy losses of b and c quarks in the region 
of hard collisions |i| > A^, where the perturbative QCD 
is no applicable in the vicinity of Landau pole \t\ ~ A^. 

In the work were found the physical parameters T 
and r for expanding plasma at ^/s = 200A GeV with ini- 
tial meanings Tq = 410.6 MeV and tq = 0.405m;;^ We 
use here the number of this parameters for calculations 
of for b and c quarks at different energies by for- 
mula ([T4|) . These calculations were fulfilled by means of 
numerically. We show in the Table U the energy losses of b 
and c quarks for bg and eg collisions for series of T and r 
for expanding plasma (for example at pb.c — 8 GeV). We 
have here at \t\rnin ~ A^ for b quark {Mb — 4.75 GeV) 
Cimin = 0.31, yniin = and for c quark (M, = 1.3 



We use 



= 0.028 



0. 00645;, 

and Umi, 



0.002_Ee 



and for Landau 



GeV) ^1 mill yj.yj^^^ o-iiva yniin Cl^ 

pole Zniin — 1. The invariant analytic charge ^ in for- 
mula ([T?)) describe also the asymptotic for large z, and 
we use the ordinary QCD constant A = 200 MeV. In the 
third column of the Table H] we give for bg collisions, 
in the fourth column we give the same energy losses for 
eg, and then we give the summary energy losses {^■)tot- 

The calculation gives for summary energy losses: for b 
quark {^ft^t = 0.112, and for c quark {^Yt^t = 0.357. 

One should determine also the energy losses of heavy 
quarks at collisions with light quarks in the plasma. We 
assume that in the initial stage at T = Tq all the quarks 
are nonequilibrated. With decrease of temperature (at 
T < Tq) the number of equilibrated quarks increase. In 
the work we estimated the part /3 of equilibrated 
quarks at expansion of plasma. At T < 233 MeV the all 



quarks become equilibrated. We estimate here the heavy 
quarks energy losses also at collisions with portion of 
nonequilibrated quarks. We assume that nascent quarks 
have the similar momentum distribution as gluons (ij. 
The matrix element squared |Mp have the factor 4/9 for 
bq and cq collisions, however there is the some exceeding 
of the number of quarks over gluons one at middle tem- 
perature. We show in the Table HI] the results of energy 
losses calculation ^ for bq and cq collisions by analogy 
with the Table H] also for pb^c = 8 GeV. 

The calculation gives for summary energy losses: for b 
quark (^)L = 0.097 and for c quark {^)^^t = 0.309. 

In this work we calculated the summary energy losses 
{%)tot of b and c quarks for various energies. How- 
ever in Sec IIVI we use the fragmentation function for 
b and c quark into heavy B and D mesons with ac- 
counting of energy losses b and c quarks. Therefore one 
should evaluate the energy losses in relative values ^ 
for various p with accounting of heavy quark masses M. 
It is convenient for this to use the values where 



AE = vV+^ - V7- 
ing of the next term ^: 




AP. We have with account- 



1 - 



1 _ ^ 



(15) 

Using the calculation of summary values {^-^)tot for 
various values p and formula (1201) . we show in Fig. [5] 
the results of numeral calculations of ^(p) for bg, eg 
and bq, cq collisions {Mb = 4.75 GeV, M^ = 1.3 GeV). 
These results we use in Sec. IIVI for calculation of nuclear 
modification factor Raa- 
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IV. CALCULATION OF NUCLEAR 
MODIFICATION FACTOR Raa FOR b AND c 
QUARKS 

The heavy quarks b and c fragmentize into heavy B 
and D mesons (look at scheme ®). In the PHENIX [l| 
and STAR [S] experiments are observed only single non- 
photonic electrons in the semileptonic decays of heavy 
mesons. These fragmentation functions were considered 
in the number works, for example [l^, [23|, [HI and 
others. The fragmentation function contain usually the 
different parameters for b and c quarks. If the energy of 
heavy quark decrease before fragmentation, its transverse 
momentum is shifted by energy loss on the value A.p{p) . 
This effect consist in replace of fragmentation function 
D{z) by effective one -D{z') [13, where: 



(16) 



Ap ■ 
P 



— . We use here the Peterson fragmentation 



where z 

function [19| . which has been used most extensively in 
applications: 



NHz{l-zf 
[(l-z)2 + eQz]2- 



(17) 



The parameter Nh is normalization, and eg is ~ 

the ratio of the effective light and heavy quark masses. 
Usually is used the traditional values ~ 0.006 and 
Cc — 0.06 [i^]. However the value £& = 0.006 is appropri- 
ate only when a leading-log (LL) calculation of spectra is 
used. When NLL (the next to leading level) calculation 
are used, smaller values of e are needed to fit the data (for 
bottom quarks cross section) . In the work [IJ] is used the 
more appropriate value ef, = 0.002 (with FONLL calcu- 
lation) to fit data. In the present work we use the frag- 
mentation function pT)) with this parameter ei, = 0.002. 
It should be noted also that parameter ec — 0.06 is ex- 
tracted from LO fits to charm production data. However 
this parameter extracted from NLO fits is Ec — 0.02 [2^ . 
This parameter Cc — 0.02 we use also in the fragmenta- 
tion function (fT7|) . The spectra of heavy mesons H after 
of fragmentation of Q quark is determined by formula: 



da I 



dz 



Dq{-^) da 
z dpi' 



(18) 



da 



where p^ = / z. The heavy quarks distributions -^-tq 

for b and c quarks are shown in Fig. [TJ If the heavy quark 
Q lose Ap{p) before fragmentation, the spectra of heavy 
mesons is determined by formula: 



da2 



dpf 



^dz Dq{z') da 



(19) 



The semi-leptonic decays of heavy mesons dominate 
the electrons spectrum at RHIC experiment up to p± ~ 

20 GeV, and we use the value z^-m — ^ for every mean- 
ings oi p^. The value Zmax is determined from condition 
z' — — Jap = 1. Using this condition and Fig. [2] 

for one can build the functions Pbips) and Pcipo) 
and to find Zmax for various meanings of pb and po- As 
we are interested in the total factor R!'^^{ph), one should 
determine the values Zmax for the same momentum pB 
and pd of heavy mesons. We write down for example the 
several values of Zmax^ 

a) for bg and eg collisions: 
Pb = 6 GeV z^ax - 0.857 
Pd^Q GeV Zmax = 0.658 
PB = 6.9 GeV z„,ax = 0.86 
PD = 6.9 GeV z,nax = 0.678 
PB = 8.5 GeV w = 0.85 
PD - 8.5 GeV z„,ax = 0.71 

b) for bq and cq collisions: 
PS = 6 GeV zmax - 0.88 
PD GeV Zinax = 0.69 
PB = 6.9 GeV z„,ax = 0.864 
PD = 6.9 GeV z„,ax = 0.701 
PB = 8.5 GeV zn^ax = 0.867 
PD = 8.5 GeV z„,ax = 0.74 

The spectra of heavy mesons H after of fragmenta- 
tion without and with energy losses of heavy quarks is 
determined by formulas (I18p and (|19|) with improving 
parameters = 0.002 and ec — 0.02 of Petersons func- 
tion. The normalization factors Nb and Nd is fixe d by 
the probability that b quark fragment into B meson [l^ : 
1 

J dzD^{z) = 1 and by analogy for . Also the nor- 

malization factors N'g and N'jj for modified fragmenta- 
tion functions is fixed by condition [221] : 



dz' 



1 



= 1 



(20) 



and similarly for D^{z'). However it should be noted 
that for example the all b quarks with energy losses at bg 
and bq collisions come from the same spectra there- 
fore one should fix the normalization by condition: 



N'b 



dz- 



(z') 



dz- 



(z') 



P Jbq 



iz) 



(21) 

We write down here the function D^{z') without nor- 
malization factor N'g, and from here we determine nor- 
malization factor N'g. From analogous condition for c 
quarks (with D^{z') and corresponding values z^^^ we 
determine the factor N'jj). These conditions correspond 
to value of the nuclear modification factor without en- 



ergy losses - R 



b+c 
AA 



1. Thus the total factor R/Hp") 
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for b and c quarks with energy losses it is convenient to 
write in the form: 



N'gjAbg + Abq)+ N'^jA^g + A,g) 



NbAi, + NdAc 
We use here the designation in numerator: 



(22) 



A, 



bg 



1 - 



and Abg correspond to Abg^bq- 
We use also: 



Ar 



dz D^{z') da 



eg 



and Acq correspond to Acg^cq- 
The designation in denominator is 



A,^ 



and 



Ac 



dz da 
—■0^(2) — . 
z ap^ 



In these formulas the fragmentation functions Df, and 
Dc we write down also without normalization factors. 

However in experiments are observed only non- 
photonic electrons in the semileptonic decays of heavy 
mesons. It was noted 0], that apparently the electrons 
from b mesons dominate the spectrum for p± > 4 GeV. 
But recently the contribution of B meson decays to non- 
photonic electrons inp+p collisions at -^=200^ GeV has 
been measured using azimuthal correlation between non- 
photonic electrons and hadrons The extracted B de- 
cay contribution is approximately 0.5 at p± > 5 GeV/c. 
The ratio is found to be 0.52 ± 0.03. We use 

this result for calculation of factor i?^)^ by formula , 
assuming the absence of the exceeding of electrons spec- 
trum from B mesons. We use the formulas (|2T|). and ((22|) 



for calculation of factors N'g and N'j^ and also of nuclear 
modification factor I^AAiP^) ■ Figures [1] 

and [5] for numeral calculation on the z integrals for ev- 
ery of fixed meanings p^ — p^. We calculate also the 
corresponding factors Nb and Nd without energy losses. 
We show here the results of numeral calculation of fac- 
tor i^^iA (PJ-)^ for example at pf^pf^ 6.9 GeV/c in 
accordance with designations in formula (j22p for every 
member in this formula: N'g = 14 55+14 147 : ^bg + A^q — 

0.172 + 0.175, iVl, = 2.646+2.742 ^ ^cg- 

and Nb ^ TJ 



2.646+2.742 ' 

Au = 0.277, Nr 



Acq = 0.02 + 0.025, 
Ac = 0.18. 



Table III: The nuclear modification factor for collision h and 
c quark with gluons only and with gluons and light quarks 



(GeV/c) 


RAA^o 


nQia+i) 


5.5 


0.22 


0.235 


6 


0.235 


0.251 


6.9 


0.265 


0.281 


8 


0.295 


0.315 


8.5 


0.315 


0.335 



3.377 



In result we have the nuclear modification factor for 
pf = 6.9 GeV/c from formula R''/a{p1 = 6.9 

GeV/c) = 0.281 

We give in the Table Hill the results calculations of fac- 
tor i?^^^ for some meaning p^ at collisions of b and c 
quarks with gluons only (i?^^) and with gluons and light 

quarks(i?^i^+^^) 

We see that accounting of 6 and c quark collisions also 
with light quarks gives exceeding of value R''^^ in less 
than 10 per cent. We show in Fig. [3] the comparison 
of nuclear modification factor i?AA (Pj^) at 5.5 GeV/c 
< p^ < 8.5 GeV/c in the present model with experimen- 
tal data for central Au + Au collisions at y/s = 200 A 
GeV [I]- We have agreement in afore-cited region oi p^. 
The not great increase of b!'^^ with increase of p^ ap- 
parently is also in experiment. We do not calculate this 
factor i?^^ for too low meanings of p^, where the Pe- 
tersons fragmentation functions in form (jl7p apparently 
inapplicable [l9l |. 



V. CONCLUSION 

We investigate in this work the collision energy losses 
of b and c quarks at hard collisions with gluons and light 
quarks in QGP produced in central Au + Au collisions at 
RHIC energy. We assume that heavy quarks not come 
to an equilibrium with partons in plasma. The hard col- 
lisions are described by formulas of perturbative QCD 
(pQCD) at large enough momenta transferred. The hard 
collisions in QCD limited by momenta ~ A^, where in 
the pQCD we have Landau pole in the coupling constant 
as{Q^) at Q2 = (the A ~ 200 MeV is dimensional 
constant of QCD). In this work we use the analytic mod- 
els for the running coupling in QCD, which satisfied to 
conditions of analyticity in local Quantum Field Theory 
(look at Sec.ini. We use here the analytic model 3 
with invariant charge in the form Q , which have no Lan- 
dau pole and describe also the asymptotic at large mo- 
menta transferred > A^. In Sec. IIIII we calculate the 
energy losses of b and c quarks at hard collisions with 
light partons at momenta transferred > A^. We show 
in Fig. [2] the results of calculations of summary losses 
Ap/p of b and c quarks with momentum p(, c at hard 
collisions with gluons and also with light quarks. These 
results we use in Sec. IIVI for calculation of the total nu- 
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clear modification factor R^a- ^^^^ fragmentation 
function for b and c quarks into B and -D mesons with ac- 
counting of summary losses Ap/p and also fragmentation 
functions without energy losses. We use here the Peter- 
sons fragmentation function [l^ , which has been widely 
adopted in analysis determining the " hardness" of heavy 
quarks fragmentation functions. Its attraction being that 
it has only one free parameter, eg, for each heavy quark 
Q. But we use the Petersons fragmentation functions 
with improving parameters eg for b and c quarks to fit 
the data, as has been mentioned in Sec. IIVI We take 
into account to the equal degree the energy losses of b 
and c quark at calculations. One should not neglect by 
b quark contribution, thought it is less than c quark one. 
We take into account also the relative contribution of B 
decays ~ 50 percents in non-photonic electrons We 
show in Fig [3] the comparison of factor R^^A ^'^^ heavy 
mesons B and D in the present model with experimental 



data [^. These data show the analogous suppression of 
nonphotonic electrons in the region of middle meanings 
of p±. It is interesting to compare the energy losses of 
heavy quarks and nuclear modification factor with oth- 
ers models. We plan also the analogous calculation with 
others fragmentation functions. 
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Figure 1: The differential cross section of charm (the upper curve) and bottom (the lower curve) quarks in pp collisions at 
= 200 GeV from [3], 



Figure 2: The relative the energy losses Ap/p b and c quarks at various momenta. The upper curves correspond to c quark for 
eg collisions(the curve 1) and for eg + cq (the curve 2). The lower curves correspond to analogous collisions for b quarks (the 
curve 3 and 4) 



Figure 3: The comparison of nuclear modification factor Rj^^ in the present model at middle meaning of p± (the solid curve) 
with experimental data of STAR collaboration ^] for nonphotonic electrons. 
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